Microscale fluid flow induced by thermoviscous expansion along a traveling wave 
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The thermal expansion of a fluid combined with a temperature-dependent viscosity introduces 
nonlinearities in the Navier-Stokes equations unrelated to the convective momentum current. The 
couplings generate the possibility for net fluid flow at the microscale controlled by external heating. 
This novel thermo-mechanical effect is investigated for a thin fluid chamber by a numerical solution 
of the Navier-Stokes equations and analytically by a perturbation expansion. A demonstration 
experiment conflrms the basic mechanism and quantitatively validates our theoretical analysis. 
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Spatial confinement of a liquid changes its fiow be- 
havior markedly, since the importance of surface forces 
relative to the volume forces increases as the confinement 
becomes smaller. Recently, flow at the scale of millime- 
ters and below has attracted signiflcant attention, stimu- 
lated by the rapid advances to manipulate and to control 
small-scale devices 0, S 3 • Since microfluidic flow of- 
ten is essentially at zero Reynolds number, viscous drag 
overwhelms the inertial effects of the fluid giving rise to 
peculiar flow behavior [3, @|- In particular, heat conduc- 
tion becomes quite efficient, since the thermal diffusiv- 
ity K = k/pcp, with k the thermal conductivity and pcp 
the volumetric heat capacity, implies thermal relaxation 
times /tt^k of the order of 0.1ms for water confined 
by walls separated hy D = 10 pm. Usually, this strong 
thermal coupling gives rise to a uniform temperature and 
all physical processes are isothermal. Since pressure ef- 
fects upon density are negligible if the velocities are small 
with respect to the speed of sound, fiuid fiow is essentially 
"dynamically incompressible" i-e. the velocity field is 
solenoidal, div v = 0. 

Recently, it has been noted by Yariv et al. [§| that for 
the case of unsteady heating this is no longer the case 
and one can in principle generate a non-solenoidal fiow. 
A transient fiuid fiow emerges due to a thermal expansion 
of the fiuid as a reaction to the heating in combination 
with the thermal diffusion. 

In this Letter we propose a novel mechanism to gen- 
erate net fiow in a thin fiuid chamber, i.e. a viscous 
liquid confined between two plates separated by a dis- 
tance of the order of a few micrometers. The driving of 
the fiuid fiow is provided by imposing a traveling tem- 
perature wave. We show analytically within a thin-film 
approximation that such unsteady heating leads to net 
fiuid fiow. Then we corroborate our analytic approach 
by a finite-element calculation. Last we provide first ex- 
perimental evidence that there is indeed net fiow and 
that the fundamental dependences have been correctly 
identified. 

The basic mechanism may be summarized as follows: 



Due to thermal expansion of the liquid the motion of 
the heating results in a pressure modulation. The pres- 
sure gradients induce a potential fiow which, however, 
does not lead to net fiuid fiow. Yet, the tempera- 
ture dependence of the shear viscosity gives rise to a 
small net mass transport typically opposite to the mo- 
tion of the heat source. The fiow velocity is then propor- 
tional to the thermal expansion coefficient of the fiuid, 
a = pQ^{dp/dT) as well as to the thermal viscosity coef- 
(Lcient (3 = -r]^\dTj/dT). 

The geometry under consideration consists of a thin 
chamber of height D and a much larger lateral exten- 
sion. There is a natural small parameter e = D/L <^ 1, 
where L is a typical length scale in the lateral direction. 
Since the film is thin one expects the Navier-Stokes equa- 
tions to be dominated only by a few terms. To iden- 
tify the relevant contributions dimensionless quantities 
are introduced as t ~ Tt*,Vi = Uv*,Vi = vv^,Xi ~ 
Lx*,x± = Dx*^,p — Pp*. Lateral directions are la- 
beled by a latin index i G {x,y}, whereas the vertical 
direction is indicated by the subscript _L. To keep all 
terms in the mass conservation law the scales are cho- 
sen as V = eU,T = L/U, and consequently the fiuid 
fiow is essentially in-plane. To balance the leading or- 
der term in the momentum conservation law the pres- 
sure scale has to be chosen as P = riQU/{t^L). Momen- 
tum transport via convective processes may be ignored 
from the very beginning since the Reynolds number is 
small, fie = poUL/r]a ^ 1 [l3|. Then the mass conser- 
vation law and the momentum balance in perpendicular 
and lateral direction read to leading order in a thin-film 
approximation, see e.g. (1, l^, ll| (restoring units) 
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As a consequence of the thin film geometry, (i) bulk vis- 
cous processes do not contribute to the leading order 
equations, (ii) the pressure is homogeneous in the per- 
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pendicular direction, (iii) lateral pressure gradients drive 
the lateral fluid flow. Note that the divergence of the 
velocity field is irrelevant for the momentum balance but 
plays a crucial role in the mass conservation law. 

The coupling to the temperature enters the equations 
via the expansion of the fluid as well as by the tem- 
perature dependence of the shear viscosity. For the 
problem at hand it is appropriate to neglect the me- 
chanical compressibility of the fluid and to expand the 
equation of state to first order in the temperature field, 
p = po(l ~ a5T), where 5T is the local temperature 
change and a ~ —{din p/dT)p the thermal expansion 
coefficient at the reference state {po,To). The density 
field then is eliminated from the equations of motion in 
favor of the temperature field ST. A second ingredient 
for the theoretical model is to include the variation of the 
shear viscosity with temperature. Introducing the ther- 
mal viscosity coefficient /3 = —{dhir]/dT)p, the coupling 
reads rj = i]o{l — (3ST) to first order. 

In the case of considerable thermal coupling to the 
walls the temperature field approximately refiects the 
profile of the heat source and then the temperature is 
uniform perpendicularly to the walls. The in-plane vari- 
ations of the temperature field occur on scales A, 6 » £>, 
where A denotes the wavelength and b the typical lateral 
extension of the heating. Consequently, the spatial de- 
pendence of the density and the shear viscosity is only 
in-plane as it is inherited from the temperature profile. 
With this additional assumption, the in-plane momen- 
tum balance equation, Eq. (fTc|) . is readily integrated for 
no-sHp boundary conditions at the walls 
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= - — z(D - z)V^p, 
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i.e. the velocity profile corresponds to Poiseuille fiow. 

Since the spatial dependence of the velocity profiles 
perpendicular to the wall is known, one may project the 
three-dimensional problem to the two-dimensional plane 
by averaging over the vertical z-direction. 
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Because the density is assumed to vary only in the lat- 
eral direction, the mass conservation law, Eq. (fTa|) . can 
be directly averaged over the perpendicular direction, 
dtp + '^iipvi) = 0, without introducing new terms. The 
pumping process becomes stationary in the frame of ref- 
erence comoving with the heat wave, i.e. substituting 
Vx ^ Vx — u, where u denotes the velocity of the heat 
wave. Then in the comoving frame the averaged mass 
conservation law yields 



uVxP + pViVi = -ViVip . 



(4) 



Here the velocities Vi are small quantities, typically Vi ^ 
u. From the closed set of equations, Eqs. ([3]) and Q, the 
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FIG. 1: (color online) Finite element solution for a tempera- 
ture wave in the comoving frame, (a) The temperature pro- 
file is indicated by the color scale and the arrows refer to the 
induced velocity flow, (b) The corresponding pressure modu- 
lation (color scale) and the solenoidal component of the flow. 



fiuid velocity can be calculated. The temperature profile 
in the comoving frame enters Eq. (jH in the form of den- 
sity gradients, Vip = —apoViST and Eq. ^ via the tem- 
perature dependent viscosity. The net flow emerges only 
if the thermal viscosity coeflicient /? is non- vanishing. For 
a temperature-independent viscosity the pressure acts as 
a velocity potential, Eq. ([3]), and averaging along the di- 
rection of the wave propagation yields a vanishing net 
flow. 

We exemplify the flow profiles for the case of a gaussian 
temperature wave (5T(7^ = STq cos{kx) , obtained 

by a numerical solution of Eqs. |[3]) and ^ using the Fem- 
lab(§)3.1 software package by Comsol. We have used a 
wavelength A = 27r/fc = 50 pm, a lateral width b = 7 pm 
of the gaussian modulation, a temperature amplitude of 
STq = 3.8 K, and a thermal wave velocity u = lOOmm/s. 
The parameters for a = 3.0 ■ IQ-'^K'^, (3 = 0.022 K"^ 
correspond to water at ambient temperature. The pa- 
rameters have been chosen such that the temperature 
modulation appears approximately as a train of circular 
heated areas mimicking the demonstration experiment 
below. The regions where the temperature gradient is 
maximal (minimal) act as sources (sinks) for the veloc- 
ity field, giving rise to a locally dipolar pattern of the 
fiuid fiow, see Fig. [T][a). Careful inspection reveals that 
the velocities in the heated regions are bigger than in the 
colder ones. The pressure profile is shifted by a quarter of 
a wavelength reflecting the local expansion of the fluid. 



3 



Fig. W[h) . The solenoidal part of the velocity field has 
been obtained as the difference of two Femlab calcula- 
tions for thermal viscosity coefficients (3^0 and /3 = 0. 
Close to the center of the wave this contribution is di- 
rected opposite to the thermal wave velocity with an aver- 
age pump velocity v^^^ = 2.88iim/s at the center of the 
thermal wave, whereas far away from the heating there is 
a characteristic backfiow. The pressures induced by the 
pumping are below 10~^ bar (Figure lb). The pressure 
limitations of the glass chamber used are estimated to 
10 bar by a finite element calculation. However, the pres- 
sure induced by the pumping increases with the viscosity 
of the fiuid and eventually will stall the pump motion 
since thermal expansion will deform only the chamber 
walls instead of triggering fiuid motion. So we expect to 
equally pump fiuids with viscosities 10^ fold higher than 
water. 

To gain further insight, we develop a perturbative an- 
alytical solution for the fiuid fiow appropriate for small 
temperature changes. Formally, the expansion is per- 
formed in the small parameters a and /?, and we shall 
show that net fiuid fiow first occurs at 0{a • (3). The 
quiescent fiuid Vi = 0,p = po solves the governing equa- 
tions to zeroth order in a, /3, i.e. if all coupHngs to the 
temperature field are ignored. To leading order in a, the 
term on the r.h.s. of Eq. ([4]) is of second order and may 
be ignored. Similarly, ([3]) is expanded to first order in 
the temperature change 
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Already at this point one concludes that to order 
the velocity corresponds to potential fiow only, hence no 
net pumping arises to that order. The leading order to 
the pump velocity is then expected to be of order 0{a(3). 
To identify the net fiuid motion it is favorable to intro- 
duce the velocity potential, x, and the stream function, -0, 
viavx = VxX + VyV'i'^y = ^yX-^x^- From Eq. ^ one 
infers that the leading term corresponds to pure gradient 
fiow, i.e. X = — (-C^/12?7o)p, where x(^) denotes the ve- 
locity potential and f~{x,y). To the required order x 
is determined by combining Eqs. ([4]) 
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Once the velocity potential is determined, the stream 
function is calculated by extracting from (3dT{r)Vix{r) 
its solenoidal component. By Helmholtz's decomposition 
theorem, the stream function satisfies a Poisson equation 
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The lowest order contribution to the stream function is 
thus proportional to q;-/3, which also sets the overall scale 
of the thermo- mechanical effect. Furthermore one infers 
from Eq. ((Tj) , that a strictly one-dimensional temperature 



modulation does not give rise to solenoidal fiow. The 
remaining task is to solve the set of Poisson equations, 
which can be easily implemented by numerical methods. 
To determine the average fiow for a travelling periodic 
temperature wave train of wavelength A, it is actually 
sufficient to know the velocity potential. 
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where the last relation follows from the representation, 
Eq. ^ and an integration by parts. 

To illustrate the physics we consider a gaussian tem- 
perature wave ST{f) = STq cos(fca;)e~^ Z^'' , where b 
characterizes the lateral width of the wave train. The 
Poisson equation, Ecl {l]), can be solved exactly in terms 
of error functions [I4I, however, it is instructive to con- 
struct an approximate solution for wide waves, kb ^ 1. 
Then the problem is effectively one-dimensional, and x 
depends on the lateral coordinate y only parametrically. 
One readily calculates x = — ua^To sin(A:x)e~^ Z^'' /k, 
implying a net fiuid motion 



(9) 



typically opposite to the motion of the traveling tem- 
perature wave. Since by assumption the changes in the 
density and the viscosity are small, the net fiuid motion 
is much slower than the velocity of the temperature wave, 

Vfl <^ u. 

We have confirmed the predicted fiuid movement in a 
demonstration experiment. A circular variant of a ther- 
mal wave is imposed with infrared light to a thin fiuid 
film, see Fig. [21 The fiuid movement is recorded us- 
ing micrometer-sized fiuorescent particles and the ther- 
mal wave is imaged stroboscopically with temperature- 
sensitive fiuorescence. Within experimental errors, the 
theory captures the thermally triggered net fiow. 

The details of the experiment are as follows: A fiber 
laser at 1455 nm and a maximal power of 5 W (RLD-5- 
1455, IPGLaser) is defiected by an acousto-optical defiec- 
tor (Pegasus Optik, AA.DTS.XY.lOO) and moderately 
focused from below (Thorlabs, C240TM-C, / = 8 mm, 
NA = 0.5) to a 10 pm thin fiuid film sandwiched between 
sapphire windows. The light is absorbed by water with 
an attenuation length of 305 ym. The chamber is imaged 
from above by a fiuorescence microscope (Zeiss, Axiotech 
Vario) and a CCD camera (SensiCam QE, PCO). The il- 
lumination is provided by a green LED (LXHL-LX5C, 
Luxeon) which for the temperature imaging was modu- 
lated with a bandwidth of 150 kHz using a laser current 
source (LD-3565, ILX Lightwave). The temperature field 
was imaged using 50 yM of the fiuorescent dye BCECF 
under stroboscopic illumination (isl . 14 1 to determine the 



temperature amplitude of the thermal wave. 
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FIG. 2: (color online) Demonstration experiment, (a) A cir- 
cular thermal wave is created in water by a moving focus of 
an infrared laser. The temperature is inferred by stroboscopic 
imaging of a temperature-sensitive fluorescent dye. (b) Peak 
temperature along the circumference. A sinusoidal fit yields 
a temperature amplitude of STo = 3.8 ± 0.5 K. 
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FIG. 3: (color online) Experiment versus theory, (a) The 
fiuid velocity scales quadratically with the amplitude of the 
thermal wave STq. (b) The speed increases linearly with the 
velocity it of the thermal wave. Measurements are given as 
dots, the theory as solid lines with the geometrical prefactor 
determined in Ref. [lal. Error bars are standard errors. 



with the thermal wave, opposite to the observed flow di- 
rection. Under the experimental conditions, both effects 
would yield similar attractive forces. However, control 
experiments in heavy water with its 100-fold decreased 
absorptive heating only revealed moderate attraction into 
the illuminated ring, but no pumping movement along it 
(<0.13 ym/s). This indicates that the particles are ideal 



tracers of the fluid motion at the used focus size 

With the above described temperature pattern and a 
thermal wave velocity of m = 100 mm/s we measure a 
peak fluid flow of = 2.7 ± 0.5pm/s. We expect 
from the wide wave solution, Eq. [9l a fluid velocity of 
^theo _ 4.77jini/s. Applying the complete analytical so- 
lution [ill or the Femlab result yields v^^^ = 2.88 ym/s 
which describes the experimental result quantitatively. 

Two parameter variations further conflrm the theoret- 
ical model (Fig. [S]). The dependence on the temperature 
amplitude STq exhibits a parabolic dependence (Fig. [3] a) 
as expected theoretically by Eq. ([9]). Furthermore, the 
fluid velocity scales Hnearly with the velocity u (Fig.[3]b) 
of the thermal wave as predicted (solid lines) . We there- 
fore flnd experimentally that a circular thermal wave trig- 
gers a fluid flow according to the theoretical description. 

To conclude, we showed that a thermal wave can move 
a fluid by the nonlinear combination of the temperature- 
dependent density and viscosity. While the obtainable 
fluid velocities are small at present, future improvements 
could allow microfluidic applications. 

We thank Joseph Egger for initial discussions. F.M.W. 
and D.B. were supported by the Emmy Noether Scholar- 
ship of the DFG. Support by the Nanosystems Initiative 
Munich (NIM) is gratefully acknowledged. 



The thermal wave is generated by a scanning laser 
spot. Six individual temperature peaks are created by 
a circular pattern of radius R= 50 ym with a base fre- 
quency of 30 kHz. The frequency is faster than the ther- 
mal relaxation time, which is dominated by the verti- 
cal heat currents and determined to 0.21ms using a fl- 
nite element calculation. As result, a steady thermal 
wave of six points is generated, measured in Fig. [51 
The points are rotated using a slower shift frequency 
= 0.2 ... 2 kHz. The result is a circular thermal wave 
with velocity u = fs x 2nR/6 = 10 . . . 100 mm/s. 

Fluid velocities are measured by single particle track- 
ing of 80 pM of 1pm diameter silica beads (PSi-Gl.O, 
Kisker). Particle velocities could be well discriminated 
from Brownian diffusion by tracking 5 independent par- 
ticles over 10pm with a positional error of 1 pm. The 
peak velocity of the parabolic flow proflle was measured 
by selecting the beads in the chamber center. The exper- 
imental values were multiplied by 2/3 to compare with 
the chamber-averaged theoretical values. Optical trap- 
ping or thermophoresis would move the tracer particles 
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